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The longitudinal asymmetry induced by parity- violating (PV) components in the nucleon-nucleon 
potential is studied in the charge-exchange reaction 3 He(n,p) 3 H at vanishing incident neutron en- 
ergies. An expression for the PV observable is derived in terms of T-matrix elements for transitions 
from the 2S+1 Lj = ^o and 3 Si states in the incoming n- 3 He channel to states with J = and 
1 in the outgoing p- 3 H channel. The T-matrix elements involving PV transitions are obtained in 
first-order perturbation theory in the hadronic weak-interaction potential, while those connecting 
states of the same parity are derived from solutions of the strong-interaction Hamiltonian with the 
hyperspherical-harmonics method. The coupled-channel nature of the scattering problem is fully 
accounted for. Results are obtained corresponding to realistic or chiral two- and three-nucleon 
strong-interaction potentials in combination with either the DDH or pionless EFT model for the 
weak-interaction potential. The asymmetries, predicted with PV pion and vector-meson coupling 
constants corresponding (essentially) to the DDH "best values" set, range from -9.44 to -2.48 in 
units of 10 -8 , depending on the input strong-interaction Hamiltonian. This large model depen- 
dence is a consequence of cancellations between long-range (pion) and short-range (vector-meson) 
contributions, and is of course sensitive to the assumed values for the PV coupling constants. 
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PACS numbers: 21.30.-x,24.80.+y,25.10.+s,25.40.Kv 



INTRODUCTION, RESULTS, AND 
CONCLUSIONS 



A number of experiments aimed at studying parity vio- 
lation in low-energy processes involving few nucleon sys- 
tems are being completed or arc in an advanced stage 
of planning at cold neutron facilities, such as the Los 
Alamos Neutron Science Center, the NIST Center for 
Neutron Research, and the Spallation Neutron Source at 
Oak Ridge. The primary objective of this program is to 
determine the fundamental parameters of hadronic weak 
interactions, in particular the strength of the long-range 
part of the parity- violating (PV) two-nucleon (TV TV) po- 
tential, mediated by one-pion exchange (OPE). While 
such a component is theoretically expected on the ba- 
sis of the weak interactions between quarks and the 
spontaneously-broken chiral symmetry of QCD, exper- 
imental evidence for its presence has proven to be elu- 
sive, and indeed current constraints are inconclusive, for 
a review see Ref . [l[ . 

In contrast, in the strong-interaction sector OPE domi- 
nates the NN potential at internucleon separations larger 
than 1.5 fm, and the spatial-spin- isospin correlations it 
induces leave their imprint on many nuclear properties. 
These include, for example, i) the observed ordering of 
levels in light nuclei and, in particular, the observed ab- 
sence of stable systems with mass number A = 8 [2J , ii) 
the single-particle energy spacings and shell structure of 
medium- and heavy-weight nuclei [3| and, in particular, 
the observed change in the energy gap between the /111/2 
and <?7/2 orbits in tin isotopes [4j, and hi) the relative 
magnitude of the momentum distributions of pp versus 



np pairs in nuclei [5j , which leads to the strong suppres- 
sion of (e, e'pp) relative to (e, e'np) knock-out cross sec- 
tions from 12 C, recently measured at Jefferson Lab Q. 

The determination of the parameters that character- 
ize parity violation in nuclei requires evaluating matrix 
elements of hadronic weak-interaction operators between 
cigenstates of the strong-interaction Hamiltonian. Thus, 
experiments in this field are especially reliant on the- 
ory for their analysis and interpretation. For this rea- 
son, over the last several years, we have embarked on 
a program aimed at developing a systematic framework 
for studying PV observables in few-nucleon systems, for 
which accurate — essentially exact — calculations are pos- 
sible. Two earlier papers [7|,|8| dealt with the two-nucleon 
system, and provided a rather complete analysis of the 
longitudinal asymmetry in p-p scattering |7| up to 300 
MeV lab energies, and of a variety of PV observables in 
the np system [8J, including, among others, the neutron 
spin rotation in ft-p scattering and the photon angular 
asymmetry in the n-p radiative capture at thermal neu- 
tron energies. In the next phase, we have studied the 
spin rotation in ft-d [9( and n-a [10( scattering at cold 
neutron energies. 

Measurements are available for the following PV ob- 
servables: the long itudinal analyzing power in p-p [llj- 
J14j and p-a [15| scattering, the photon asymmetry 
and photon circular polarization in, respectively, the 
^(n.i^H [Il]-[l3 and ^(n^H pj] radiative cap- 
tures, and the neutron spin rotation in ft-a scatter- 
ing J13 . [2fJ. There is also a set of experiments which 
are currently being planned, including measurements 
of the neutron spin rotation in n-p [l9( and n-d [2l| 



scattering, and of the longitudinal asymmetry in the 
charge-exchange reaction 3 He(n, p) 3 H at cold neutron en- 
ergies [221, tnc subject of the present paper. 

At vanishing neutron energies, the only channels en- 
tering the incoming n- 3 He scattering state have quan- 
tum numbers 2S+1 L,j = x So and 3 Si. In the out- 
going p- 3 K scattering state, the relevant channels are: 
2S+1 Lj = 1 So, 3 Si, 3 Di with positive parity, and 3 P , 
1 Pi, 3 Pi with negative parity. We show (in Scc.|TT| that 
the PV observable in this process, i.e. the longitudinal 



analyzing power A z , reads 
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where 6 is the angle between the proton momentum and 
the neutron beam direction, and the coefficient a z can 
be expressed in terms of products of T-matrix elements 
involving (three) parity-conserving (PC) and (three) PV 
transitions as 
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In T L g L , s , the label J specifies the total angular mo- 
mentum, the superscripts 21 denote the charge-exchange 
transition n- 3 He to p- 3 H (as opposed, for example, to 
the elastic transition, which would be denoted by the su- 
perscripts 22), the subscripts LS (L'S') are the relative 
orbital angular momentum and channel spin of the n- 
3 He (p- 3 H) clusters, and lastly the overline is to note the 
inclusion of a convenient phase factor — see Eq. ([To]) be- 
low. The PC (PV) T-matrix elements have L + L' even 
(odd), and the sum E in Eq. ([3]) is proportional to the 
3 Hc(n,p) 3 H cross section. We observe that a z vanishes if 
only the channels 1 So and 3 Po (with J = 0) are retained. 
The T-matrix elements arc related to the (real) R- 
matrix elements (Sec. IIIII and Appendix [A"]) , and the 
latter for PC transitions are calculated via the Kohn 
variational principle with the hypcrsphcrical-harmonics 
(HH) method [23l,|24} (Sec.EJ- We use strong-interaction 
Hamiltonian models, consisting of the Argonnc vis 
(AV18) Hl| or chiral (N3LO) [H] two-nucleon potential 
in combination with the Urbana IX (UIX) [27| or chiral 
(N2LO) |28[ three-nucleon potential. The HH calcula- 
tion is a challenging one, for two reasons. The first is the 
coupled-channel nature of the scattering problem: even 
at vanishing energies for the incident neutron, the elas- 
tic n- 3 He and charge- exchange p- 3 H channels are both 
open. The second is the presence of a J" = + reso- 
nant state (of zero total isospin) between the p- 3 H and 
n- 3 He thresholds, which slows down the convergence of 
the expansion, and requires a large number of HH basis 
functions in order to achieve numerically stable results. 
Further discussion of this aspect of the calculations is in 
Sec. |V1 where we also present current predictions for the 
n- 3 Hc scattering lengths corresponding to the Hamilto- 



nian models mentioned earlier. They are in good agree- 
ment with the measured values. 

The i?-matrix elements involving PV transitions are 
computed in first-order perturbation theory with Quan- 
tum Monte Carlo techniques (Sec.lVip. We adopt as PV 
potential the meson-exchange (DDH) model of Desplan- 
ques et al. [29|] as well as the pionless effective- field-theory 
(EFT) model recently derived in Refs. (MIMl (Scc. HV|) , 
and present results for the various components of the 
DDH and EFT potentials in combination with the AV18, 
AV18/UIX, N3LO, and N3LO/N2LO Hamiltonians in 
Sec. I VIII Additional results for the R- and T-matrix 
elements, and combinations thereof entering the PV ob- 
servable, are listed (for the AV18/UIX) in Appendix [B] 
for completeness. For the DDH model only, we also 
present predictions for a z corresponding essentially — but 
see Sec. IIVI for further details — to the "best values" of 
the 7T-, p-, and oj-meson weak-interaction coupling con- 
stants [29(. These predictions range from -9.44 to -2.48 
in units of 10 -8 depending on whether the N3LO/N2LO 
or AV18/UIX Hamiltonian is considered, and thus ex- 
hibit a significant model dependence due to cancellations 
(or lack thereof) between the pion and vector-meson con- 
tributions. 



It is useful to express the asymmetry as 
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where the /i„'s, a = n, p, ui and i = 0, 1, 2, denote 
the PV coupling constants in the DDH model along with 
the isospin content of the corresponding interaction. The 
coefficients C* are listed in Table [I] and depend on the 
input Hamiltonian used to generate the continuum wave 
functions, as well as on the assumed values for the PC 
pion- and vector-meson coupling constants and associ- 
ated cutoffs (see Table HV| . 



J 



AV18 

AV18/UIX 

N3LO 

N3LO/N2LO 



-0.1892(86) 
-0.1853(150) 
-0.1989(87) 
-0.1110(75) 



C» 
-0.0364(40) 
-0.0380(70) 
-0.0120(49) 
+0.0379(56) 



+0.0193(9) 
+0.0230(18) 
+0.0242 (9) 
+0.0194 (10) 



CI 
-0.0006(1) 
-0.0011(1) 
(-0.0002(1) 
-0.0007(1) 



-0.0334(29) 
-0.0231(56) 
+0.0080(30) 
+0.0457(36) 



Cl 
-0.0413(10) 
-0.0500(20) 
-0.0587(11) 
-0.0408(14) 



TABLE I: The coefficients C* entering the PV observable a z , corresponding to the AV18, AV18/UIX, N3LO, and N3LO/N2LO 
strong-interaction Hamiltonians. The statistical errors due to the Monte Carlo integrations are indicated in parentheses, and 
correspond to a sample consisting of ~ 130k configurations. 
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The coefficients G™, follow from the linear combination 
given in Eq. flSJ). Isotensor p-exchangc (C7 2 ) is negligible. 
The isoscalar and isovector vector-meson exchanges give 
contributions of the same magnitude, both of which are 
smaller than OPE. However, the OPE contribution seems 
to be significantly suppressed. For example, in the case 
of the neutron spin rotation in n-d scattering this contri- 
bution is calculated to be at least a factor of ~ 30 larger 
than that of any of the p and ui exchanges, which is not 
the case for the process under consideration. This may 
be due to the predominant isoscalar character of the '■So 
and 3 Po channels — see discussion in Appendix [Bj The 
N3LO/N2LO results should be considered as preliminary, 
since the HH solution for the + wave function has not 
yet fully converged (at least as far as the singlet scat- 
tering length is concerned, sec Sec. fV|. This fact may 
explain wh y th e inclusion of a threc-nucleon potential 
like N2LO |28[ should reduce C\ by almost a factor of 
two relative to the other models. This point will be dis- 
cussed in Sees. IV1 and IVII Finally we note that the "best 
values" for the PV couplings constants of the pion and 
p-meson are (in units of 10 -7 ) respectively +4.56 and 
-16.4, and this leads to the large cancellation (and con- 
sequent model dependence) in the values predicted for a z 
and referred to earlier. 

We conclude by observing that the EFT analysis pre- 
sented in this work could be improved by employing chi- 
ral potentials in both the strong- and weak-interaction 
sectors. At order Q/A x , where Q is the low en- 
ergy/momentum scale that characterizes the particular 
process of interest, and A x ~ 1 GeV is the chiral- 
symmetry-breaking scale, the PV potential contains 7 
low-energy constants (LECs), 5 of which are associated 
with four-nucleon contact terms, and the remaining 2 
with long-range OPE components [3(|. When electro- 
magnetic interactions are also introduced, another (un- 
known) LEC must be included — it is needed to fix the 
strength of a PV two-body current operator of pion 
range [30(. One can envisage, at least in principle, a 
suite of experiments involving A = 2-5 systems, which 
would constrain, in fact over-constrain, these eight LECs. 
Some of these have been mentioned above, additional 
ones include, for example, measurements of the photon 
asymmetries in the radiative captures 2 H(n, 7) 3 H and 



3 Hc(n, 7) 4 He. These processes are strongly suppressed: 
the experimental values for the corresponding (PC) cross 
sections [32|, [33( are, respectively, almost 3 and 4 orders 
of magnitude smaller than measured in 1 H(n, 7) 2 H. One 
would naively expect relatively large PV asymmetries in 
these cases, possibly orders of magnitude larger than in 
the A=2 system. Clearly, accurate theoretical estimates 
for them could be useful in motivating our experimen- 
tal colleagues to carry out these extremely challenging 
measurements. 

From a theoretical perspective, most of the method- 
ological and technical developments needed to carry out 
the calculations are already in place. We have recently 
reported results [3J] for the A = 3 and 4 (PC) captures, 
using wave functions obtained from the N3LO/N2LO 
Hamiltonian and electromagnetic currents derived in chi- 
ral EFT up to one loop [35[ , which are in excellent agree- 
ment with data. However, there is one aspect in the com- 
putation of the proposed PV threshold captures, which 
still needs to be addressed: the determination of the 
small admixtures induced by the PV potential into the 
bound and continuum wave functions. Even a first-order 
perturbative treatment of those admixtures requires con- 
struction of the full Green's function for the strong (PC) 
Hamiltonian, an impractical task. However, it may be 
possible to generate them using correlated basis meth- 
ods, similar to those employed in Ref. (36| . 



II. THE PARITY- VIOLATING OBSERVABLE 

The neutron energies in the reaction 3 Hc(n, p) 3 H of in- 
terest here arc in the mcV range, and at these energies 
only two channels are open: the n- 3 He elastic channel and 
the p- 3 H charge-exchange channel. In the following, the 
index 7=1 (2) is used to identify the p- 3 H (n- 3 He) clus- 
ters in the final (initial) state. In the absence of strong 
and weak interactions between the two clusters, the wave 
function in channel 7 is written as 
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where 'J™ 3 is the (antisymmetrizcd) trinucleon bound- 
state wave function in spin projection m^, x^ 1 is the nu- 
cleon spin-isospin state with spin and isospin projections 
m\ and p for 7=1 or n for 7=2, respectively, and <p is the 
inter-cluster wave function, i.e. a Coulomb wave func- 
tion for 7 = 1 or simply a plane wave e lq2 ' yp for 7 = 2. 
The separation between the center-of-mass positions of 
the two clusters is denoted by y p with y p = r; — Hijk, 
and their relative momentum is specified by q 7 , so that 
the energy E is given by 
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where ^[(q^; y p ) reduces to a regular Coulomb function 
Fl (97; Up) (multiplied by a phase factor we need not spec- 
ify here) for 7 — 1 or a spherical Bcsscl function x jl (x) 
for 7 = 2, with x = q 7 y p - 
arc defined as 
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while the Clebsch-Gordan coefficients associated with the 
re-coupling of the angular momenta (and other factors) 
are lumped into 



Here B 1 and M 1 are the binding energy and mass of 3 H 
( 3 He) for 7=1 (2), and m 1 is the proton (neutron) mass 
for 7=1 (2). Lastly, the wave functions in Eq. ([5]) arc 
antisymmetrizcd by summing over the four permutations 
p with (ijk,l) = (123,4), (124,3), (134,2), and (234,1). 
It is useful to expand the wave functions in Eq. |JSJ| in 
partial waves as 
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The momentum q 7 has been taken to define the spin- 
quantization axis, i.e. the z-axis. 

In the presence of inter-cluster interactions, the 7i- 3 He 
wave function in the asymptotic region reads 
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and contains outgoing spherical waves in the n- 3 He elas- 
tic channel (7 = 2) as well as in the p- 3 K charge-exchange 
channel (7 = 1) multiplied by corresponding T- matrix el- 



ements T 



77', j 



lsl'S'- Here 771 = a/ii/qi, where a is the fine 
structure constant and \i\ is the p- 3 H reduced mass de- 
fined above, and ol is the Coulomb phase-shift. Thus 
Coulomb distortion in the p- 3 H outgoing state is fully 
accounted for. 

The probability amplitude M m i m i t m3 mi to observe a 
p- 3 H final state with spin projections m! x and to 3 , respec- 
tively, is obtained from 
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where we have assumed that the p- 3 H state is in partition 
(123,4) corresponding to permutation p = 1, namely the 
bound cluster consists of particles 123 and the proton is 
particle 4. For brevity, we have also set y = y p =i. Using 
the orthonormality of the channel functions f^^g „, we 
find 
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where the Clebsch-Gordan coefficients require J z = S z 
TO3 + mi , <S^ = m 3 + mi , and M' = J z — S' z = to 3 + mi 
(m 3 + to'tJ. 



The spin-avcraged cross section follows from 
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since (l/4)(gi//xi) |M m ^ mi , m3 , mi | 2 dO is the flux of 
outgoing particles in the solid angle dfi = dy, and 
(1/4) ((72/^2) is the incident flux, where the factors 1/4 
originate from the normalization factors 1/a/4 in Eqs. © 
and (fTTj) . These cancel out in Eq. (fT3|) , leaving an extra 
1/4 coming from the average over the initial polariza- 
tions. The longitudinal asymmetry A z is defined as 
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TABLE II: The coefficients Ci 1 
channels. 
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At meV energies it suffices to keep only L = in the 
entrance channel, so that 
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where the P| (<?)'s are associated Legcndre functions, 9 
is the angle of the outgoing proton momentum relative 
to the direction of the incident beam, the C/ 1 ^ 2 s (|M|)'s 
denote combinations of Clebsch-Gordan coefficients, de- 
fined as 
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After inserting the expression for Z , ' , " and carry- 
ing out the sums over mi, 7713 and 7711,7713, we find the 
unpolarized cross section to be given by 
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where in the second line we have ignored T-matrix el- 
ements involving transitions to odd parity final states 
(and hence parity violating), since these are induced by 
hadronic weak interactions and consequently are much 
smaller than the parity-conserving T-matrices associated 
with strong interactions. We observe that the matrix ele- 

ments T 21,J (and T ' ) are finite in the limit q^ = 0, and 
therefore oyj is divergent as 92 goes to zero, as expected 
for a neutron capture reaction. 

The asymmetry A z can be written as 
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and lastly the phase factor tL x L 2 
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ensures that either L\ or L2 must be odd, which in turn 
implies that either T 0J ' 2 5 or T 0J j? s involves a parity- 



violating transition, 



a transition from an incoming 



,(17) positive parity n- 3 He state to an outgoing negative parity 
p- 3 H state. The non-vanishing C's for the relevant chan- 
nels are listed in Table [U and evaluation of the sums in 
Eq. (Tl8|) allows one to express the parity-violating asym- 
metry as in Eqs.([T])-([3]). 



III. T-MATRIX ELEMENTS 

The calculation proceeds in two steps: we first deter- 
mine, via the Kohn variational principle, the i?-matrix 
elements, and then relate these to the T-matrix elements. 
The wave function describing a scattering state with to- 
tal angular momentum JJ Z in channel ~/LS is written 
as 



7,-Lo 7,-ko 



-* 



FJJ Z 
f,LS ' 



E 

7'L'S' 



DTYV yG,JJ z 

n LS,L'S'^i',L'S' 



(21) 



,\,JJ Z 



where the asymptotic wave functions ^ ' ' ' LS Z with A 
F, G are defined as 
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and the superscript A = F is to denote the regular ra- 
dial functions introduced earlier in Eq. ([?])■ and A = G 
is to denote the irregular Coulomb or spherical Bcsscl 
functions, namely 

7 = 1: T?(x) = G L ( m ,x) ; 7 = 2: ^ l g (.t) = -ifr L (s) . 

(23) 
The tilde over Gl and j/l indicates that they have been 
multiplied by short-range cutoffs in order to remove the 
singularity at the origin. Thus J^ is well-behaved in all 
space. The normalization factor Z? 7 , 
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(24) 



and k = v/3/2, is introduced for convenience — k is a 
numerical factor relating the inter-cluster separation y p 
to the Jacobi variable Xi p , i.e. xi p = ny p (see Eq. (|44]) 
below). 

The wave functions V ' LS 2 vanish in the asymptotic 
region, and describe the dynamics of the interacting 
nucleons when they are close to each other, while the 



^ls'l'S'' s are t ne -R-matrix elements. The latter, as well 
as the coefficients entering the expansion of Vt 7 ' LS Z in 
terms of hypcrsphcrical-harmonics functions, are deter- 
mined via the Kohn variational principle 



R 



77 ,J 
LS,L'S' 



W7,J 
l L'S',LS 



m^r.-(^ J r.\H-E\^ JJ 



1,LS\ 



j',L'S' 



(25) 



as discussed in Sec. [V] 

The next step consists in relating the R- to the T- 
matrix elements. To this end, it is convenient to simplify 
the notation by dropping the superscripts JJ Z and by 
introducing a single label a to denote the channel quan- 
tum numbers LS, so that the wave functions in Eq. (|2ip 
corresponding to 7 = 1 and 2 are written as 



\T> =<I»f + vJ/;f 



■E^*w+EC f v -( 27 ) 



From these we form the linear combination 

* = Yl (Vet,*'®!,*' + V a , a >$2, a >) 



(28) 



where the matrices U and V arc determined below. In- 
serting the expressions above for ^7,0 and rearranging 
terms lead to 



J 



* 



* ( 



n 



Y\U-i (UR 11 + VR^)] a Jl a , + Y (UR 11 + VR^) aQ , (*&, + i *[ a , 



V - i {UR 11 + VR 22 ) 



* 



a, a' 



E (^ 12 



vi? 



22\ 



* 



im 



(29) 



r 



where ^ c is a combination of internal parts of no interest 
here. We now require ^ to consist, in the asymptotic 
region, of a plane wave in channel 7=2 (or n- 3 Hc) and of 
a purely outgoing wave in channel 7=1 (or p- 3 H). These 
requirements are satisfied by demanding that 



U - i {UR 11 
V - i {UR 12 



VR 21 ) = 



(30) 
(31) 



VR ZZ ) = I , 
where / is the identity matrix. Comparing the resulting 



^> with the wave function given in Eq. (jT0|) . specifically 
its component in channel LSJ, allows one to express the 
T-matrix as 



T 



21, j 

LS,L'S' 



Jh_ 

D 2 q\ 



{U J R 



7 nil,,/ 



J r>21,./\ 



V J R 



1 LS.L'S' 



D 2 qi 



U 



LS.L'S 1 ' 



(32) 



where we have reinstated the LSJ labels. Finally the 
matrix U is obtained by solving the system in Eq. pip : 



J 



r 21.J 
1 LS.L'S' 



D 2 qi 



I - i R 22 ' J + R 21 ' J (I - i i? 11 - 7 )" 1 ^ 12 " 7 

I 



R^il-iR 11 ^)- 1 



J LS,L'S' 



(33) 



In fact, we compute the i?-matrix elements at zero en- ergy, i.e. in the limit q 2 — > 0, and define 



pl2,J 
-pl2,J _ n LS.L'S' 

^LS,L'S' ~ 17+572 ' 

92 



■7521, J 



R 



R 



LS,L'S' 



21, J 
LS,L'S' 
£+1/2 ' 
92 



10 

11 

12 



+ 



gTr h^ 
2^/2 m 
9ph°„ 



9pfe"(l + «p) 

m 

9p h P 

2m 

9p fc p( 1 + K p) 

2 m 

9ph 2 



2V6m 

3 P fe p( 1 + «p) 
2\/6m 



9wh£(l+Ku) 



2 771 

flu) ftl(l+K„) 



2m 



/r H (o 

A(r) 
/p(r) 
A(r) 

/pW 
fp(r) 
fp(T) 

Mr) 
Mr) 
Mr) 

Mr) 
Mr) 
Mr) 



2 ji" 



C 6 



o 



£(C 2 + C 4 ) 

X 



2 M 2 r 



2)i' 



2)i' 






/n EFT (r) 
/,W 




Mr) 


Mr) 



/,W 








o 



Tjrr 



(Ti X Tj) z ((Ti + (Tj) X 
T; ■ Tj (o-j - O-j) ■ X (2) 
Tj • Tj (fTi X (Tj) ■ X 



(1) 
*Ji- 



(3) 
»Ji- 



(Ti + Tj) z {(Ti - (Tj) ■ X$ + 

,(5) 



(-r i + r J ),(o- l x ov,-) ■ X^ ; _ 



\0 Ti,z"Tj,z 
\0 Ti^z Tj , z 



n ■ T, 



(<Ti 



Ti ■ Tj) ((Ti X (Tj) ■ X 



v(6) 
(7) 



((Ti — (Tj 

{(Ti X O-j) -X 



,(8) 
Sj,4 

(9) 

»J>- 



X 



(10) 



,(11) 



(Ti+Tj)4^ XO"j)-X^_ 
(12) 



(TJ - Tj)« (<Ti + ffj) • XJ 



TABLE III: Components of the DDH and EFT models for the parity-violating potential. The vector operators X'™_ and 
functions fx(r), x = n, p, to, fx, are defined in Eqs. (|38[) - (|39|) and Eqs. (|40[) - (|41[) . respectively. Only 5 operators and low-energy 
constants enter the pionless EFT interaction at the leading order, and in this paper they have been chosen to correspond to 
the rows 1, 4, 6, 8 and 9. 



R 



■22, J 
LS,L'S' 



R 



22. J 

LS.L'S' 



(34) 



and it can be shown that the i?-matrix elements are finite 
in this limit. In particular, we note that the factor q\ fol- 
lows from the small argument expansion of the spherical 



Bessel function jl in V ' 2 LSi 



1/2 



LS , while the extra q 2 is due 
to the normalization Z?2- At zero energy, we have 



I-iR 22 - J + R 21 ' J (I-tR ll ' J )- 1 R 



12, J 



/ 



(35) 



One is the model developed thirty years ago by Desplan- 
ques et al. |29| (and known as DDH): it is parametrized 
in terms of 7T-, p-, and w-meson exchanges, and involves 
in practice six weak pion and vector-meson coupling con- 
stants to the nucleon [37|. These were estimated within 
a quark model approach incorporating symmetry argu- 
ments and current algebra requirements [29l . |38| . Due to 
the inherent limitations of such an analysis, however, the 
coupling constants determined in this way have rather 
wide ranges of allowed values. 



since R 22,J and the product R 21 ' J R 12 ' J are proportional 
to g2 or higher powers of qi . Furthermore, the relevant T- 
matrix elements entering the expression for the asymme- 
try A z are those with quantum number L = in channel 
7 = 2, and hence 



T, 



21, J 



OJ.L'S' 



s/oi hi,, 



21, J 
0J,L"S" 



[I-iR 



11. .A 



'L"S".L'S' 



L"S" 



(36) 

with J — 0,1. Note that we have neglected the difference 
in the n- 3 He and p- 3 H reduced masses. 



IV. THE PARITY- VIOLATING POTENTIAL 

Two different models of the PV weak-interaction po- 
tentials are adopted in the calculations reported below. 



The other model for the PV potential considered in 
the present work is that formulated by Zhu et at I30| in 
2005, and reduced to its minimal form by Girlanda [31| 
in 2008, within an effective-ficld-theory (EFT) approach 
in which only nucleon degrees of freedom are retained 
explicitly At lowest order Q/A x , where Q is the small 
momentum scale characterizing the low-energy PV pro- 
cess and A x ~ 1 GeV is the scale of chiral symmetry 
breaking, it is parametrized by a set of five contact four- 
nucleon terms. 



J 



The DDH and EFT PV two-nucleon potentials are con- 
veniently written as 



E c «° 



(n) 



a = DDH or EFT 



(37) 



where the parameters c" and operators O 



in) 



12, are listed in Table HTT1 In this table the vector 



1, 

operators X^"V are defined as 



ij,± 



r{n) 






Pij ; Jn{ r ij)\ + 
Pij 5 Jnyij )\_ 



(38) 
(39) 



where [...,.. 
commutator (- 
ator, p^ = (pj- 
f x (r), x = TT,p 



_. denotes the commutator (— ) or anti- 

) , and pij is the relative momentum oper- 

-Pj)/2. In the DDH model, the functions 

and u, are Yukawa functions, suitably 



modified by the inclusion of monopole form factors, 



/*(r) 



1 

Ait r 



A x r 
2 



™ x 
A 2 

1 1 ■., 



(40) 
In the EFT model, however, the short-distance behavior 
is described by a single function f p (r), which is itself 
taken as a Yukawa function with mass parameter fi, 



47rr 



(41) 



with p ~ TOtt as appropriate in the present formulation, 
in which pion degrees of freedom are integrated out. 

In the potential wB DH , the strong- interaction cou- 
pling constants of the n-, p-, and w-meson to the nu- 
cleon are denoted as g^, g p , k p , g u , Ku, while the weak- 
interaction ones as h\, h° pl h„, h 2 pl h%, h^, where the su- 
perscripts 0, 1, and 2 specify the isoscalar, isovector, and 
isotensor content of the corresponding interaction compo- 
nents. In the EFT model, the five low-energy constants 
C\, C\, C2 + C4, C5 and Cq completely characterize v F J FT , 
to lowest order Q/A x . 



ffa/47T 



13.2 

0.840 
20.0 



6.1 
0.0 



W 



-16.4 
3.23 



10' x hi 



4.56 

-2.77 
1.94 



10' x hi 



-13.7 



A a (GeV/c) 



1.72 
1.31 
1.50 



TABLE IV: Values used for the strong- and weak-interaction 
coupling constants and short-range cutoff parameters of the 
7T-, p-, and o;-meson in the DDH potential. 

The values for the coupling constants and short-range 
cutoffs in the DDH model arc listed in Table llVl while the 
mass p in the EFT model is taken to be m^. These values 
for coupling constants and cutoffs were also used in the 
DDH-based calculations of PV two-nucleon observables 
in Refs. |7|, |8[ and neutron spin rotation in rid scatter- 
ing [9j. In particular, we note that the linear combination 



of p- and w-meson weak coupling constants correspond- 
ing to pp states has been taken from an earlier analysis 
of pp elastic scattering experiments |7| . The remaining 
coup ling s are the "best value" estimates, suggested in 

Ref. daf. 

In the analysis of the a z observable to follow, we will 
report results for the coefficients Z^ DH and /^ft in the 
expansion 



12 



E« ■ 



(42) 



Thus we will not need to consider specific values (or range 
of values) for the strength parameters c". However, the 
I" depend on the masses (and short-range cutoffs A x 
for the DDH model) occurring in the Yukawa functions. 
Note that the coefficients G % a entering Eq. (|4|) are ob- 
tained from the J° DH 's and c° DH 's listed in Tabic |HT] 



via 



cl 



C° p 



cl 



cl 



>h 



-ii 



DDH 



2^21 

9p ,ddh 
m 



9 P (1 + k p ) L 

m 

gp(l + K p ) 

2m 
9p0- 



DDH 

3 



DDH 

5 



Cl = 



2m 4 

9p jDDH 

2\/6m 6 2y/6i 

9u r DDH 9uQ- + Ku) r DDH 
"~ i 8 ~ 1 9 



2m 



[ 12 



-I L 7 



DDH 



(43) 



9u rDDH _ 9 u ( 1 + k oj ) ,DDH _ am tL) 

2 to 10 2 m n 2 m 12 



THE HH WAVE FUNCTIONS 



The "internal" wave function *£'/s% see Eq. (|2T]). is 
expanded in the HH basis. For four equal mass particles, 
a suitable choice for the Jacobi vectors is 



xi p 

x 2p 
x 3p 




(44) 



where p specifics a given permutation corresponding to 
the ordering (ijkl). By definition, the permutation p = 1 
is chosen to correspond to (1234). 

For the given Jacobi vectors, the hyperspherical coor- 
dinates include the so-called hyperradius p, defined by 



- . U2 



>lp 



b 2p 



3p 



(independent of p) 



(45) 



and a set of angular variables which in the Zernike and 
Brinkman [39j, |40j representation are (i) the polar angles 
Xi p = (9i P ,(j>ip) of each Jacobi vector, and (ii) the two 



additional "hypcrspherical" angles (f>2 P and </>3 P , defined where Xj P is the magnitude of the Jacobi vector x Jp . The 

set of angular variables Xi p ,X2p,X3 P ,02 P , and </>3 P is de- 
noted hereafter as J7 p . A generic HH function reads 



as 

COS 4> 2 p 



X 2p 



,2 
'lp 



''2;j 



COS( 



':!;; 



%3p 



'lp 



C 2p 



r 2 
! 3/; 

(46) 



H 



KAM 



(fi„) 



JX, n 2 n 3 



Y £l (xi p ) ® Yi 2 (x 2p ) ® ^3 (x 3p ) 
£2 



(sin (/> 2p ) f J (cos </> 2p ) 2 (sin </> 3p ) f 1_ 



2+2«2 



Ail/ 



x (cos ^ P£ +1/2 ■ £2+1/2 (cos 20 2p )P^+ 2 " 2+2 ■ ^ +1 /2 (cos 20 3p ) , 

I 



(47) 



where P£' b are Jacobi polynomials, and the coefficients 
A/"„2„ 3 are normalization factors. The quantity K = £±+ 
^2+^3 + 2 («2 + "3) is the so-called grand angular quan- 
tum number. The HH functions are the eigenfunctions 
of the hyperangular part of the kinetic energy operator. 
Another important property is that p K ^^ 2 ^ L2 „ 2 „ 3 (ftp) 
are homogeneous polynomials of the particle coordinates 
of degree K. 

A set of antisymmetrizcd hypcrangular-spin-isospin 
states of grand angular quantum number K , total orbital 



angular momentum A, total spin S, and total isospin T 
(for the given values of total angular momentum J and 
parity n) can be constructed as follows: 



* 



KAT.T 



XX AET 

p=l 



(ijkl) 



(48) 



where the sum is over the 12 even permutations p = ijkl, 
and 



$ 



if AST 



(ijkl) 



n 



KAM 

3-^2^12^3 



(ftp)® 



J 



lx»®Xjlo ®Xk 



s,, 



®Xi 



-/■/: 



[6®&] T ®& 



(49) 



TT = 



Here, x« (CO denotes the spin (isospin) state of particle i. 
The total orbital angular momentum A of the HH func- 
tion is coupled to the total spin S to give the total angu- 
lar momentum JJ Z , whereas the parity it is (— X)*- 1 2 3 ■ 
The quantum number T specifies the total isospin of the 
state, and /i labels the possible choices of hyperangular, 
spin and isospin quantum numbers, namely 



H = {£i,£2,is,L2,n2,n,3,Sa,Sb,T a ,Tb} 



(50) 



compatible with the given values of K, A, S, T, J, 
and 7r. Another important classification scheme for the 
states is to group them in "channels": states belong- 
ing to the same channel have the same values of angular 
(^1,^2,^3,^2, A), spin (S a ,Sb,T,), and isospin (T a ,T b ,T) 
quantum numbers, but different values of n 2 and 713. 

Each state -$)KAY,t en tering the expansion of the four- 
nucleon wave function must be antisymmetric under the 
exchange of any pair of particles. Consequently, it is 
necessary to consider states such that 



<I> 



ifAST 



(ijkl) 



-$ 



if AST 



(jlkl) 



(51) 



which is fulfilled when the condition 

£3 + S a +T a = odd , 



(52) 



is satisfied. 

The number Mkat.t of antisymmetrizcd functions 
^ , ^ AET having given values of K, A, S, and T, but differ- 
ent combinations of quantum numbers \x — see Eq.(|50j) — 
is in general very large. In addition to the degeneracy 
of the HH basis, the four spins (isospins) can be coupled 
in different ways to total £ (T). However, many of the 
states \D f ™ ET , with /1 ranging from 1 to Mkajyt, ar e 
linearly dependent. In the expansion of V " LS Z , it is nec- 
essary to include only the subset of linearly independent 
states, whose number is fortunately significantly smaller 
than Mkat,t- 

The internal part of the wave function can be finally 
written as 



* 



C,JJ Z 

l-LS 



E 

ifAST 



E 



-y.LS 
ifASTp, 



(p)K 



if AST 



(53) 



where the sum is restricted only to the linearly indepen- 
dent states. We have found it convenient to expand the 



10 



"hypcrradial" functions u^ AET (p) in a complete set of 
functions, namely 



M-l 



II 



■y,LS 



-KAT.T^P) = Yl C MET ftm 9m{p) , (54) 

and have chosen 



m=0 



>0>) 



•(8) 



111 



(m + 8)! 



^ 2 ^(Me- W2 



(55) 



The c coefficients of the expansion (f5"4"| and the R- 
matrix elements of Eq. pip are determined variationally 
via the Kohn variational principle. This principle states 



where L\ ' (ftp) are Laguerre polynomials |4l| 



that the functional 



nil ,J 



defined in Eq. (|25|) is sta 



tionary with respect to variations in the RJJs'l'S' ano - 



1,LS 



By applying this principle, 



a linear set of 



-KA£Tn,m- 

equations for RJ/s'l'S' ano - c mstu m ^ s obtained [23| . 
then solved using the Lanczos algorithm. The other pa- 
rameter entering the expansion is the (non linear) pa- 
rameter j3 (see Eq. (|55|) ). used to describe the hyperra- 



dial functions 



1,LS 



'■ka'ETu(p)- We have checked that, once 
a sufficient number M of functions g m (p) are employed 
(M « 20), the results are practically independent on j3. 
In the present work we have used /3 = 4 fm . 

The application of the method has two main difficul- 
ties. The first is the accurate computation of the matrix 
elements of the Hamiltonian. By exploiting the proper- 
ties of the HH functions, however, this task can be notice- 
ably simplified, as discussed in Refs. [2j,|421- The second 
difficulty is the slow convergence of the HH expansion. 
This problem has been overcome by dividing the set of 
states iJ/^" AET defined in Eq. (j48f (in the following re- 
ferred to simply as "HH states") in classes, depending 
on the value of C = l\ + £2 + £3, total isospin T, and 77,2 
and 713. In the present paper, we have considered four 
different classes. Since for n- 3 He scattering the asymp- 
totic states do not have a definite total isospin (they are 
a superposition of T = and T = 1 components), it is 
mandatory to include HH states with both T — and 1. 
The contribution of T = 2 states is expected to be tiny 
and consequently they have been ignored in the present 
paper. 

Following Refs. [42, 53[i in the first class we have in- 
cluded the n 2 = HH states belonging to some special 
channels, for which the convergence has been found to 
be critical. The radial part of these HH states depends 



only on cos < 



'3p 



j/p and thus they take into account 



two-body correlations. The n 2 > HH states belonging 
to the same channels are included in the second class, 
together with those having C < 2. The other classes are 
then defined simply by grouping HH states belonging to 
channels with an increasing value of C. In particular, for 



the construction of the positive (negative) parity "inter- 

C 1 7 7 

nal"wave function ^ ' L g z , classes 3 and 4 include all HH 
states with C = 4 and 6 (£ = 3 and 5), respectively. 
The convergence of these last two classes is less critical, 
and consequently, only HH states with lower values of 
grand angular quantum number K need be considered. 
Moreover, the convergence with C is quite fast. In par- 
ticular, we have found that, at the energy considered, the 
contribution of HH states with C > 6 can be neglected. 

The calculation is performed including in the expan- 
sion all HH states belonging to classes i = 1,...,4 
with grand angular quantum number K < Ki, where 
K\, . . . ,K<± are a set of nonnegative integers. The con- 
vergence of a quantity of interest (for example, the phase- 
shifts, or the coefficient a z defining the PV asymmetry) 
is then studied by increasing the values of Ki. A more 
complete study of the convergence will be presented else- 
where [Hj]. 

To exhibit the convergence pattern, we report in Ta- 
ble |V] the calculated n- 3 He scattering lengths. As is evi- 



(56) 



dent from Eq. (fTOjl . they are defined as 



aj 



lim T, 



22, J 



0,7,0.7 ' 



with both incoming and outgoing 7i- 3 He clusters in 
relative S-wave. Note that in general this scattering 
length is complex, since the channel p- 3 H is always 
open, and therefore the unitarity condition imposes that 
Imaj < 0, since the total cross section is proportional to 
J2j=o 1(2^ + 1) I m ^oj'oj- T ne resm ts obtained for the 
singlet (J = 0) and triplet (J = 1) scattering lengths are 
reported in Table [V] for all four potential models used 
in this work. The calculated n- 3 He scattering lengths 
are compared with experimental values and the results 
of other calculations available in the literature. 

Inspection of the table shows that the convergence for 
the triplet scattering length is very good, and that there 
is reasonable agreement with available experimental val- 
ues, and the results of other calculations, in particular 
those of the AGS method. In the case of the singlet 
scattering length, the situation is more delicate, since in 
the channel J v = + the n- 3 He interaction is attrac- 
tive and the wave function must be orthogonal to the 
4 Hc bound state. Consequently, the convergence is more 
problematic, in particular for the N3LO/N2LO interac- 
tion model. In the row labeled "EXT" , we have reported 
the extrapolated values for this quantity obtained by an- 
alyzing the convergence pattern. For the AV18, N3LO, 
and AV18/UIX interaction models we observe reasonable 
agreement with the results of other calculations and the 
experimental data. The N3LO/N2LO values are signifi- 
cantly different from those obtained with the other inter- 
action models, which is presumably related to the slow 
convergence observed in this case. A complete study of 
the n- 3 He scattering lengths is in progress |24| . 



J 



11 



Triplet scattering length <n (fm) 


Ki 


K 2 


K 3 


Ki 


AV18 


N3LO 


AV18/UIX 


N3LO/N2LO 


28 
30 


28 
30 


20 
22 


20 
22 


3.56- i 0.0078 
3.56 -i 0.0077 


3.47 -i 0.0047 
3.46 -i 0.0048 


3.39 -i 0.0059 
3.39 -i 0.0059 


3.37- i 0.0042 
3.37- i 0.0042 


RGM [44] 
FY [45] 

AGS [46] 


3.45 -i 0.0066 
3.43 -i 0.0082 
3.51 -i 0.0074 


3.56 -t 0.0070 
3.47 -i 0.0068 


3.31 - j 0.0051 
3.23 -i 0.0054 





R-matrix 
EXP [47] 
EXP [48] 
EXP [49] 



441 



3.29 

3.28(5) 
3.36(1) 

3.48(2) 



■i 0.0012 
■i 0.001(2) 



Singlet scattering length ao (fm) 


Ki 


K 3 


K 3 


Ki 


AV18 


N3LO 


AV18/UIX 


N3LO/N2LO 


48 
50 


44 
46 


30 
32 


22 

24 


7.34 -i 6.27 
7.41 -i 6. 16 


7.38- i 5.23 
7.40- i 5.20 


7.90 -i 3.65 
7.90 -i 3.59 


4.45 -i 9.02 
5.25 -i 9.25 


EXT 


7.69 -i 5.70 


7.57- i 4.97 


7.89 -i 3.44 


6.02 -i 9.48 


RGM [44] 
AGS [46] 


7.78 -i 5.02 
7.80 -i 4.97 


7.82- i 4.51 


7.62 -i 4.09 





R-matrix [44] 
EXP [47] 
EXP [48] 
EXP [49] 



7.40 - i 4.449 

7.37(6) -i 4.448(5) 

7.46(2) 

7.57(3) 



TABLE V: Convergence of the n- He singlet and triplet scattering lengths corresponding to the inclusion, in the internal part of 
the wave function, of four different classes in which the HH basis has been subdivided. For the singlet scattering length, the line 
labeled "EXT" reports the extrapolated values obtained by examining the convergence pattern with increasing number of HH 
functions in the expansion. The calculated scattering lengths are compared with results obtained using the Resonating Group 
Method (RGM), Faddeev-Yakubovsky (FY) equations, Alt-Grassberger-Sandhas (AGS) equations, as well as with results of 
R-matrix analyses. The experimental values are reported in the rows labeled "EXP" (the imaginary parts are taken from 
Ref. El). 



r 



Recently, there has been a new measurement [50| for 
the quantity Aa' = Re(ai — ao) = —4.20(3) fm. The cal- 
culated values of Aa' with the AV18, N3LO, AV18/UIX, 
and N3LO/N2LO models are -4.13, -4.11, -4.50, and 
—2.65 fm, respectively. Again the N3LO/N2LO value 
stands out: it is off that obtained with the other interac- 
tion models and the measured value. 

The convergence for the negative-parity states is sim- 
ilar to that discussed above. For the 0~ state, there is 
a close resonant state and the convergence is slow as in 
the + case. For the 1~ state, the resonance is far and 
we observe good convergence, as for the 1 + state. Note, 
however, that in these cases the N3LO/N2LO conver- 
gence pattern is not different from that observed with 
the other models. 



VI. CALCULATION 

There is a total of two (four) states with J = 
(J = 1): one (two) with positive parity having LS = 00 
(LS = 01, 21) and one (two) with negative parity having 

LS = 11 (LS = 10, 11). The .R-matrix elements R^s'^s 
with LS = 00 or LS = 11 for J = 0, and RjJsvs' with 



LS, US' = 01, 21 or LS, L'S' = 10, 11 for J = 1, involv- 
ing parity-conserving transitions induced by the strong 
interactions are calculated with the HH method, as de- 
scribed in the previous section. However, the .R-matrix 
elements involving parity-violating (PV) transitions are 
obtained in first-order perturbation theory as 



R 



7l',J 
LS,L'S' 



(* 



JJz 

7', L'S' 



PV 



* 



JJ r \ 
1,LSI 



(57) 



where L + L' must be odd. Specifically, the i?-matrix 
elements relevant for the calculation of the asymmetry 

are: R 0Q \i and i? 0011 for J = 0, and i? i'iO' ^oi'iij 



R 



ii,i 



,ii,i 



21.10' -"-21,11 



R n 



21 1 

l0 , and -Rqi'ii f° r J = 1- Quantum 
Monte Carlo (QMC) techniques are employed to evaluate 
these matrix elements (see below). 

The asymmetry in Eq. ([2]) is expressed in terms of T- 
matrix elements, which are in turn derived from .R-matrix 
elements via Eq. (|36[) . This latter equation can be further 
simplified by retaining only linear terms in the PV R- 
matrix elements, and the resulting expressions for the 
PC ^oo,bo' r oi,'oi an d ?oi,2i an d PV T Q0 ' n , T 01 \ Q and 
T 01 \ x matrix elements are listed in Appendix [A"l 

The QMC techniques used to evaluate the matrix ele- 
ment in Eq. (J5T|) are similar to those discussed in Ref. Q 
for the neutron spin rotation in rid scattering. The wave 
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functions for an assigned spatial configuration specified 
by the set of Jacobi variables (xi , x 2 , X3 ) are expanded on 
a basis of 16 x 6 spin-isospin states for the four nuclcons 

as 



ij<; 



-0(xi,x 2 ,x 3 ) = 2jVa(xi,x 2 ,x 3 ) I a) 



(58) 



where the components i/) a (xi,x 2 ,X3) are generally com- 
plex functions, and the basis states | a)= \ (n i)i(p 4- 
h(nl)3{pl)4,), \(nl)i(nl)2(pl)a(pl)4), and so on. Ma- 
trix elements of the PV potential components are written 
schematically as 



(f\0\i) = ^2 / dxidx 2 dx3^ a (xi,X2,X3) 

a, 6=1"' 

X [0(xi,X2,X3)] a6 Vi,6(Xl,X2,X 3 ) , (59) 



where [0(xi,x 2 ,x 3 )] a& denotes the matrix representing 
in configuration space any of the components in Table Hill 
Note that the operators X 2 -"' occurring in vf v are con- 
veniently expressed as 



r(n) 



-i[2/„(ry)V i3 -+?y/;(r y )] , (60) 



X i7 ._ — fjj fnVij) J 



(61) 



where the gradient operator Vjj = (V, — V_,)/2 acts on 
the right (initial) wave function, and f'(x) = df(x)/dx. 
Gradients arc discrctized as 

Vi, Q -0(xi,x 2 ,X3) ~ [ip(. ..ri+Se a ...) 

-i>(...n-6e a ...)]/(2 6),(62) 

where d is a small increment and e a is a unit vector in the 
a-direction. Matrix multiplications in the spin-isospin 
space are performed exactly with the techniques devel- 
oped in Ref. [5l[. The problem is then reduced to the 
evaluation of the spatial integrals, which is efficiently car- 
ried out by a combination of MC and standard quadra- 
tures techniques. We write 



(/I O \i)= / dxidx 2 dx3F(xi,x 2 ,x 3 )~ — 



N c 



F(c) 



(63) 
where the c's denote collectively (uniformly sampled) di- 
rections xi and Jacobi coordinates (x 2 , x 3 ), and the prob- 
ability density W(c) =| *(x 2 ,x 3 ) | 2 /(Air)— *(x 2 ,x 3 ) is 
the triton bound-state wave function normalized to one — 
is sampled via the Metropolis algorithm. For each such 
configuration c (total number N c ), the function F is ob- 
tained by Gaussian integrations over the x\ variable, i.e. 



F{c) 



E 

a,b=l 



dxi xfipf Jx 1 ,x 2 ,x 3 ) 



x [O(xi,x 2 ,x 3 )] a& -0 J ,&(xi,x 2 ,x 3 ) . (64) 
Convergence in the x integrations requires of the order 
of 50 Gaussian points, distributed over a non-uniform 
grid extending beyond 20 fm, while iV c of the order of 
a hundred thousand is sufficient to reduce the statistical 
errors in the MC integration on the PV T-matrix ele- 
ments at the few percent level. In this respect, we note 
that these errors are computed directly, by accumulat- 
ing, in the course of the random walk, values — and their 
squares — for the appropriate linear combinations of R- 
matrix elements, as given in Eqs. (|A5|) and (|A14|) (|A15|) 
of Appendix [A"l Because of correlations, the errors on the 
T-matrix elements obtained in this way are much smaller 
than those that would be inferred from the i?-matrix el- 
ements by naive error propagation. 

The present method turns out to be computationally 
intensive, particularly because of the large number of 
wave functions (and their derivatives) that have to be 
generated at each configuration (x 1; x 2 ,x 3 ). The com- 
puter codes have been successfully tested by carrying out 
a calculation based on Gaussian wave functions for the 
initial and final states, as described in the following sub- 
section. 



A. Test calculation 

In order to test the computer programs based on QMC 
techniques, we carried out a preliminary calculation using 
wave functions for which it is possible to evaluate the 
matrix elements of the PV potential also analytically. 
These (antisymmetric) wave functions are written as 



* 



JJ Z 

1,LS 



-t 



e ^ p 2 y L+2n p 



p=l 



Y L (y P )® <j>-y(ijk)®X'y(l) 



, (65) 



JJ; 



where ^> 7 (x) represents a three-nucleon (single-nucleon) 
spin-isospin one-half state with isospin projection -1/2 
(+1/2) for 7 = 1 (p- 3 H channel) and +1/2 (-1/2) for 
7 = 2 (n- 3 He channel). Thus, as in the realistic case, 
the wave functions above do not have a definite total 
isospin T but, rather, are combinations of T = and 
T = 1 states (having, of course, T z = 0). The whole 
radial dependence is given by the factor y p nfi e -f i P ; 
where p is the hyperradius. The non-negative integer np 
and the real parameter j3 can be varied so as to obtain a 
family of wave functions. For the purpose of computing 
matrix elements of two-body operators, it is convenient 
to express the pieces in Eq. (|65)) corresponding to per- 
mutations p / 1 in terms of quantities relative to the 
permutation p = 1 or (123,4). This can be accomplished 
by making use of the properties of Wigner coefficients: 
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h 1 *3 



P 'PP 2 V r L5J r" 1 

[[£i®G»] Til ®$ 



yj 8 (*3)® [xi®Xa], 



r 3 



^ 2 (X 2 )«)X3 



<g> 



5^i (xi) 8X4 



I 



(66) 



where \i an d & are the spin and isospin states of nu- 
cleon i, Xj are the Jacobi vectors corresponding to the 
permutation p — 1 and n\ + ri2 + U3 = L + 2 np , and 
the C's denote combinations of Wigncr coefficients. It is 
now relatively simple to evaluate the matrix of the PV 
potential Yli<j v Tj V = 6uf^ |9|, by expressing the wave 
functions as in Eq. ()66l) . 

As an example, we report here the results obtained for 
two J = wave functions. In Table fVll we list the values 
of the quantum numbers LSJ, and parameters np and 
/3, used in the actual calculation. The ket I7) with 7 = 
1 (2) describes a "p- 3 H" ("n- 3 He") state. We compute 
the matrix elements in two ways: i) by performing the 
analytical calculation via the transformation of Eq. (f6"fj)) , 
and ii) by using the QMC techniques discussed earlier. 



State 


r 


LS 


np 


13 


|1) 

|2) 


0+ 

0" 


00 
11 






0.25 
0.25 



TABLE VI: Values of the quantum numbers and parameters 
for some of the test wave functions used in this work. See 
text for explanation. 



The values for the matrix elements — (1|Oj!j |2) cor- 
responding to the 12 components of the DDH potential 
(see Table IIII|) are reported in Table IVIII There is good 



n 


Analytical 


QMC 


1 


-2.987 


-3.020(15) 


2 


0.000 


0.000 


3 


-0.333 


-0.349(4) 


4 


-0.264 


-0.281(4) 


5 


-0.222 


-0.233(2) 


6 


0.000 


0.000 


7 


0.000 


0.000 


8 


-0.335 


-0.349(6) 


9 


-0.143 


-0.147(2) 


10 


-0.335 


-0.349(6) 


11 


-0.286 


-0.294(3) 


12 


-0.264 


-0.281(4) 



TABLE VII: Results for the real part of the (adimensional) 
matrix element — (1| O^ |2) calculated analytically and by us- 
ing the QMC code. For the latter calculation, the statistical 
uncertainties are reported in parentheses, and correspond to 
a (rather modest) set of 5k samples. The operators O]^ are 
those of the DDH potential, listed in Table Mil 



agreement between the results of the two calculations. 
Note that the 11 = 2 contribution associated with an 
isoscalar operator as well as the n — 6, 7 contributions 
corresponding to isotensor operators vanish. The test 
wave functions consist of a superposition of T = and 
T = 1 components, and therefore it is not immediately 
apparent why this should be so. The reason for this result 
becomes clear only after carrying out the decomposition 
of the wave functions as in Eq. (|6"6"|) . It comes about 
because of delicate cancellations among various terms. 
We find it reassuring that these same matrix elements 
are seen to vanish (within machine precision) with the 
QMC code. We have verified explicitly that the close 
agreement between the two calculations persists for the 
matrix elements involving other pairs of states, including 
those having J = 1. 



VII. FURTHER RESULTS 

The results for the coefficients I" in Eq. (|4"2"|) . ob- 
tained with the (zero energy) n- 3 Hc continuum wave 
functions corresponding to the AVI 8, AV18/UIX, N3LO, 
and N3LO/N2LO strong-interaction Hamiltonians, are 
reported for the DDH and pionless EFT PV potentials 
in Tables I Villi and IIXI respectively. The subscript n in 
I" specifies the operators as listed in Table ITTTl and the 
set of cutoff parameters entering the modified Yukawa 
functions arc given in Table HVl 



A quick glance at Table IVIIII makes it clear that i) the 
contribution of the long-range component of the DDH po- 
tential due to pion exchange is at least a factor f 5 larger 
than that of any of the short-range components induced 
by vector-meson exchanges, and ii) among the vector- 
meson exchange contributions the isoscalar (n = 2, 3 and 
n = 8, 9) and isovector (n = 4, 5 and n = 8-12) ones are 
comparable in magnitude and much larger than those 
due to isotensor p- meson exchanges (n = 6, 7). It is also 
clear that the pion-exchange contribution is fairly insen- 
sitive to the choice of input strong-interaction Hamilto- 
nian (with or without the inclusion of a three-nucleon 
potential) , used to generate the n- 3 He even and odd par- 
ity states with J = and f . However, the N3LO/N2LO 
model stands out: the pion-range contribution is (in mag- 
nitude) substantially smaller than that calculated for 
the other models. Moreover, the isoscalar p-meson (u;- 
meson) contribution corresponding to n = 3 (n = 8) 
has opposite sign than obtained for the other (AV18 and 
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rDDH 


n 


AV18 


AV18/UIX 


N3LO 


N3LO/N2LO 


1 


-0.186E+00 


-0.189E+00 


-0.203E+00 


-0.113E+00 


2 


-0.826E-02 


-0.577E-02 


-0.608E-02 


-0.622E-02 


3 


+0.811E-02 


+0.864E-02 


+0.333E-02 


-0.693E-02 


4 


-0.620E-02 


-0.794E-02 


-0.970E-02 


-0.753E-02 


5 


-0.800E-02 


-0.976E-02 


-0.102E-01 


-0.781E-02 


6 


-0.359E-03 


-0.170E-03 


-0.942E-03 


+0.322E-03 


7 


+0.631E-03 


+0.115E-02 


-0.641E-04 


+0.703E-03 


8 


+0.605E-02 


+0.404E-02 


-0.699E-03 


-0.794E-02 


9 


+0.314E-02 


+0.289E-02 


-0.171E-02 


-0.577E-02 


10 


-0.689E-02 


-0.887E-02 


-0.115E-01 


-0.902E-02 


11 


-0.930E-02 


-0.113E-01 


-0.123E-01 


-0.940E-02 


12 


-0.801E-02 


-0.979E-02 


-0.115E-01 


-0.606E-02 



AV18/UIX) models. 



TABLE VIII: The coefficient 7° DH corresponding to the DDH 
potential components O^ in combination with the AV18, 
AV18/UIX, N3LO, N3LO/N2LO strong interaction Hamilto- 
nians. The statistical Monte Carlo errors are not shown, but 
are at the most 10% for the smallest contributions, and less 
than 2% for the largest. The I? DH are in units of fin -1 . 





tKFT 


1) 


AV18/UIX 


N3LO/N2LO 


1 

4 
6 

8 
9 


-0.195E+00 

-0.606E+00 

-0.639E-02 

+0.608E+00 

+0.301E+00 


-0.119E+00 
-0.391E+00 
+0.179E-01 
-0.515E-01 
+0.426E-01 



TABLE IX: The coefficient I^ FT corresponding to the pion- 
less EFT potential components O n ' in combination with the 
AV18/UIX and N3LO/N2LO strong interaction Hamiltoni- 
ans. Note that there are no potential components with n=2, 
3, 5, 7, 10, 11, and 12. The statistical Monte Carlo errors are 



not shown, but are typically less than 5%. 
units of fm~ . 



The/* 



To investigate the stability of the AV18/UIX and 
N3LO/N2LO results with respect to convergence in the 
internal part of the wave function, we present in Ta- 
ble [X] the coefficients G™ entering the PV observable a z 
in Eq. (H|) for two different choices of wave functions. 
The results labeled "wf2" were listed earlier in Table [H 
except that those relative to the N3LO/N2LO model 
are based here on a smaller number of configurations. 
These results are obtained by including in the expansion 
of the internal parts of the ± and 1 wave functions 
the maximum number of HH functions we have consid- 
ered in the present work. The results corresponding to 
the row "wfl" are obtained by reducing this number: 
in practice, for each of the classes K\ , . . . , K4 we set 
Jfi(wfl) = Ki(wi2) - 2 (see discussion in Sec. E}. Note 
also that the Monte Carlo calculation of the "wfl" co- 
efficients for the AV18/UIX model uses a factor three 
smaller number of configurations, and therefore the asso- 
ciated statistical errors are substantially larger. On the 
other hand, the "wfl" and "wf2" N3LO/N2LO results 
correspond to the same number of configurations and in- 
deed the same random walk. Taking into account errors, 
we conclude that both AV18/UIX and N3LO/N2LO cal- 
culations have converged. This is not the case as far as 
the N3LO/N2LO singlet scattering length is concerned. 





G„- 


C P U 


cj, 


c'i 


CJJ 


ct, 


AV18/UIX-wfl 

AV18/UIX-wf2 


-0.2077(281) 
-0.1853(150) 


-0.0433(116) 
-0.0380(70) 


+0.0242(29) 
+0.0230(18) 


-0.0011(2) 
-0.0011(1) 


-0.0232(77) 
-0.0231(56) 


+0.0490(30) 
+0.0500(20) 


N3LO/N2LO-wfl 
N3LO/N2LO-wf2 


-0.1118(29) 
-0.1050(35) 


+0.0369(25) 
+0.0445(33) 


+0.0200(8) 
+0.0189(9) 


-0.0009(1) 
-0.0008(1) 


+0.0390(23) 
+0.0454(31) 


+0.0402(12) 
+0.0417(12) 



TABLE X: The coefficients Ct, entering the PV observable a z , corresponding to the AV18/UIX and N3LO/N2LO strong- 
interaction Hamiltonians for two sets of wave functions (see text for details). The statistical errors due to the Monte Carlo 
integrations are indicated in parentheses. 



r 



Therefore, the differences found between the 
N3LO/N2LO and the other models are presum- 
ably due to the fact that the HH expansion for the 
N3LO/N2LO wave functions (specifically the + wave 
function) has not fully converged. Consequently, in 
the following we restrict our discussion to the results 



obtained with the AV18, N3LO, and AV18/UIX models. 
In reference to the pion contribution, the calculated 
C\ is rather insensitive to the choice of strong Hamil- 
tonian. However, there is still a considerable model 
dependence in the results obtained for the individual 
contributions due to vector-meson exchanges. This 
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model dependence, in turn, impacts very significantly 
predictions for the PV asymmetry a z , as it can be 
surmised from Table IXII Of course, this is so under 
the assumption that the values for the strong- and 
weak-interaction coupling constants characterizing the 
DDH potential arc those listed in Table HVl For example, 
the combination of coupling constants corresponding to 
pion-exchange (n = 1) and isoscalar p-meson exchange 
(n = 2 and 3) are, respectively, cf DH = (4.48 x 10~ 7 ) fm, 



„DDH 



(11.2 x 10" 7 ) fm and cf UB = (79.5 x 1CT 7 ) fm 



note that dp DH = (1 + k p )c® b}i and k p = 6.1 is the value 



adopted here for the tensor coupling of the p-meson 
to the nucleon [52j. Consequently, the contribution 
c° DH x /;P DH is comparable in magnitude and opposite 
in sign to the pion-exchange contribution cj 31311 x /p DH . 
In this respect, we note that the asymmetry a z changes 
roughly from -27 x 10~ 8 to +13 x 1(D 8 as the six PV 
weak coupling constants entering the DDH model are 
varied over their respective allowed ranges determined 
in Rcf. [23| . Thus, a z could potentially be large enough 
to make its measurement (relatively) easy. 







10 8 x 


a^ 




n 


AV18 


AV18/UIX 


N3LO 


N3LO/N2LO 


1 


-8.33±0.35 


-8.45+0.69 


-9.07+0.40 


-5.06+0.34 


2 


-9.26±0.35 


-9.09+0.70 


-9.75+0.40 


-5.76+0.38 


3 


-2.80±0.68 


-2.22+1.34 


-7.10+0.89 


-11.3+0.98 


4 


-2.86±0.68 


-2.30+1.34 


-7.20+0.89 


-11.3+0.98 


5 


-3.40±0.68 


-2.95+1.34 


-7.88+0.89 


-11.9+0.98 


6 


-3.41±0.68 


-2.95+1.34 


-7.90+0.89 


-11.9+0.98 


7 


-3.32±0.68 


-2.80+1.34 


-7.91+0.89 


-11.8+0.98 


8 


-3.97±0.69 


-3.23+1.35 


-7.83+0.90 


-10.9+0.99 


9 


-4.31±0.69 


-3.55+1.35 


-7.65+0.90 


-10.3+0.99 


10 


-4.09±0.69 


-3.26+1.35 


-7.28+0.90 


-10.0+0.99 


11 


-3.79±0.69 


-2.89+1.35 


-6.88+0.90 


-9.70+0.99 


12 


-3.45±0.69 


-2.48+1.35 


-6.40+0.90 


-9.44+0.99 


variation 


from -27.1 to +13.3 


from -27.6 to +13.8 


from -26.0 to +3.68 


from -29.7 to +6.66 



TABLE XI: Cumulative contributions to a z and associated errors (rows 1-12), obtained for the DDH PV potential with values 
for the coupling constants as listed in Table HVl The four columns correspond to the different combinations of strong- interaction 
Hamiltonians adopted in the calculations. The last row shows the minimum and maximum (central) values that a z can attain, 
as the PV couplings are varied over the allowed ranges in the original DDH formulation [291 ]. 
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The coefficients i^ FT for the operators entering the pi- 
onless EFT PV potential, that is n=l, 4, 6, 8, and 9, are 
reported in Table ITXl The coefficients I„ FT for n=l, 4, 
8, and 9, corresponding to isoscalar and isovector struc- 
tures, are all of the same order of magnitude, while that 
for n=6 with isotensor character is much smaller. Note 
that the radial functions are taken to be the same for all 
n, fn FT (r) = U(r). Of course, the I FFT, s will depend 
significantly on the value of the mass /x — cither /i = m w , 
as appropriate in the present pionless EFT formulation, 
or /j, = 1 GcV, the scale of chiral symmetry breaking, 
as appropriate in the formulation in which pion degrees 
of freedom are explicitly retained. Indeed, in this latter 
formulation the leading order component of w PV has the 
same form as the pion-exchange term in DDH. 

Finally, rough estimates have been made for the range 
of values allowed for the low-energy constants C\ , Ci + 
C4, C5, C*i, and C 6 in Ref. [3(|. However, at the present 
time a systematic program for their determination is yet 
to be carried out. In view of this, we refrain here from 



making EFT-bascd predictions for the longitudinal asym- 
metry. 
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Appendix A: From R- to T-matrices 

Consider the case with J = first. For the parity- 
conserving (PC) T-matrix we have: 



T, 



21,0 
00,00 



1 



~R 21,0 IT i pH.O'i-l 

-"-00, 00 ^ '-" Joo.oo 



11,0\-1 



+- R oo,n( / ~ i - R ' )ll 



(Al) 



where J - ji? 11 is a 2 x 2 matrix with very small off- 
diagonal elements, i.e. 



I-iR 



11,0 



a e 

e b 



— 1 R 11 ' 

, a — l i -n.00,00 ' 



• nll.O ; i • nll.O 
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(A2) 
with |a|, |6| 3> |e|. To first order in e, we approximate 

(I-iR 11 '°)- 1 = 

and hence 



1/a —e/ab \ 
-e/ab 1/6 J 



(A3) 



^21,0 _ 
J 00,00 — 



1 R 21 '° 
1 ^00,00 



qi a 



(A4) 



Similarly, for the parity-violating (PV) T-matrix element 
we find: 



rp21,0 

i oo,n 



—21,0 
00,00 



R 



(I -iR 1 



j\-i 

^00,11 



-f;21,0 



7? • IT i P n -0\-l 
■"-00,11 1 J ' n )11A1 



;7? 21 ' R 11 ' 

'-"00.00 -"-00,11 

a b 



Tj21,0 

-"-00,11 



(A5) 



The case J = 1 is somewhat more involved since the 
matrices are now 4x4. The matrix (/ — zi? 11 - 1 



=£1,1 



IS 



written as 



I-iR 11 ' 1 = 



A e 
e T B 



where A, e, and B are 2x2 matrices, 



A = 



B = 



-iR. 



i R ' 
' -"01,01 

11,1 



-iR 



01,21 
11,1 



21,01 



I-iR 

-*- R lll 



11,1 
10,10 

11, 

10 



1-iR 



-i R 11 ' 1 

' -"-10.11 

-iR 11 ' 1 

' -"11,11 



-i R 111 -i R 11 ' 1 

'-"-01,10 '-"01,11 



■iR 



21,10 



■iR 



(A6) 

(A7) 
(A8) 
(A9) 



21,11 



Note that A and B, as well as their inverse A l and B x , 
are symmetric. To first order in e, it follows that 



{I-iR 



n,i\-i 



A- 1 C 
C T B- 1 



(A10) 



where the 2x2 matrix C and its transpose are defined 
as 



C = -A- 1 eB~ 1 



C T = -B- 1 e T A' 



(All) 



This shows that (I — i R 11 ' 1 ) 1 is also symmetric in this 
approximation. The PC T 01 ' 01 and T 01 ' 21 and PV T Q1 ' 10 

and T 01 ' n matrix elements entering Eq. Q36p are then 
given by 



T, 



01,21 



+ - R 0l|2l(^ )21,01 
^Ol!oi(^~ )oi,21 



"•Eoi^lG^ )21,21 



Tn 



x(-B )l0,10 + i? i!ii(-^ )n> 



T 



21,1 
01,11 



X (B )io,n +-Rqi'ii(-B )ii,h 



(A12) 



(A13) 



-"01,01 C '01,10 + -"01,21 ( -'21,10 + -"01,10 



(A14) 



-Roi!oiCoi,n + -Roi'2lC21,ll + R 01 



10 



(A15) 



Appendix B: Numerical values for R- and T-matrix 
elements 



The set of Tables IXnUXV] arc all relative to the 
AV18/UIX+DDH model, and present results for the R- 
matrix elements involving PV transitions between states 
with J = and J = 1, the corresponding T-matrix ele- 
ments which follow from them and the parity-conserving 
(PC) R- matrix elements via Eqs. (|A5l) and (|A14|) - (|A15|) . 
and lastly the coefficients d 



(«) 



v(") 



d\ n) = T 01 ' w (n)T 



An) 



An) 



—21,0 



00,00 i 
^21,1* 



dT' = 



T 00,ll( n ) -* 01,21 
7^21,1 , ,=,21,1* 
J 01,lU n J 1 01,21 



(„) _=21,0 f 1T 21>1* 
x 2 — J 00,11 l"J J 01,01 

l 4 — J oi, ill"; J oi.oi 



(Bl) 



T 



1 



01,01 



R 



01,01 



OOoi 



where the T-matrix elements arc defined as in Eq. (|Tt 
and the label (n) on those involving PV transitions refers 



17 



n 


pll,0 

-"00,11 


-"•00,11 


1 


+0.198E+01 


+0.278E+01 


2 


+0.126E+00 


+0.305E+00 


3 


-0.149E+00 


-0.373E+00 


4 


+0.533E-01 


+0.530E-01 


5 


+0.632E-01 


+0.691E-01 


6 


+0.156E-02 


+0.154E-02 


7 


+0.129E-02 


+0.163E-02 


8 


-0.211E+00 


-0.523E+00 


9 


-0.797E-01 


-0.203E+00 


10 


+0.589E-01 


+0.588E-01 


11 


+0.720E-01 


+0.784E-01 


12 


+0.154E-01 


+0.134E-01 



TABLE XII: The parity-violating i?-matrix elements for J = 
corresponding to the DDH potential components O n ' in 
combination with the AV18/UIX strong interaction potentials 
at vanishing n- 3 He energy. The statistical Monte Carlo errors 
are not shown, but are typically ~ 1-2 % for the largest values, 
and less than 10% for the smallest. The i?-matrix element 
without (with) overline is in units of fm -1 (fm"' 2 ), see text 
for explanation. 



to the operator component O^ in Table Hill The J n 's 
discussed earlier follow from 



-1 Re 



Vzd\ 



(n) 



An) 



V2d ( J' 



+V6ct ) 



V3d { J l) 



(B2) 



where £ has been defined in Eq. <|3j> . A few words 
on units: since the operators 0^ n > do not include the 
c„'s, i.e. the combinations of nucleon mass and strong- 
and weak-interaction coupling constants, the resulting R- 
matrix (T-matrix) elements involving PV transitions are 
in units of fm -1 (adimcnsional) — they would otherwise 
be adimcnsional (in units of fm). Further, because of 
the definition in Eq. (|34|) . the R- matrix elements have 
dimensions of fm -1 / 2 . Note, however, that the T- and T- 
matrix elements only differ by a phase factor, and hence 
the former are also adimensional. 



J 



n 


R n,i 
-"01,10 


R ii,i 
-"01,11 


R ii,i 

-"21.10 


pii.i 

"21,11 


-"01,10 


-"01.11 


1 


-0.160E-01 


-0.930E-01 


+0.199E-02 


+0.365E-02 


-0.535E-01 


-0.106E-01 


2 


+0.614E-03 


+0.131E-02 


-0.216E-04 


-0.627E-04 


-0.339E-02 


+0.854E-02 


3 


-0.837E-03 


-0.198E-02 


-0.204E-04 


+0.941E-04 


+0.528E-02 


+0.604E-03 


4 


-0.188E-03 


+0.782E-03 


-0.643E-04 


-0.958E-05 


-0.179E-02 


+0.161E-02 


5 


-0.317E-03 


+0.918E-03 


-0.853E-04 


-0.133E-04 


-0.232E-02 


+0.186E-02 


6 


+0.116E-03 


+0.159E-02 


-0.404E-04 


-0.777E-04 


-0.257E-03 


-0.427E-02 


7 


-0.186E-04 


+0.191E-02 


-0.713E-04 


-0.870E-04 


+0.617E-04 


-0.518E-02 


8 


-0.769E-03 


-0.181E-02 


-0.217E-04 


+0.860E-04 


+0.506E-02 


+0.812E-03 


9 


-0.364E-03 


-0.852E-03 


-0.274E-04 


+0.382E-04 


+0.260E-02 


+0.448E-02 


10 


-0.211E-03 


+0.867E-03 


-0.716E-04 


-0.107E-04 


-0.201E-02 


+0.179E-02 


11 


-0.367E-03 


+0.105E-02 


-0.985E-04 


-0.151E-04 


-0.270E-02 


+0.214E-02 


12 


-0.543E-03 


-0.144E-02 


-0.699E-05 


+0.636E-04 


-0.258E-02 


-0.102E-03 



TABLE XIII: Same as in Table EH but for J = 1 



r 



The .R-matrix elements in J = 1 states (Table IXIII[) arc 
typically two orders of magnitude smaller than those in 
J = states (Tabic IxTl)) . Among the former, those with 
orbital angular momentum L = 2 in channel p- 3 H (7 = 1) 
arc much suppressed at the low energies of interest in 
the present work. Inspection of Table I XIII also shows 
that the (isovector) pion-exchange interaction (n = 1) is 
dominant, which suggests that the J ff = + and 0~ states 



In order to compute the d,'s in Table IXVl one needs, 
in addition to the T-matrix elements listed in Ta- 
ble IXIV| also the T-matrix elements associated with 
PC transitions. These have been calculated to be (at 
zero n- 3 Re energy): T^;° -- ■- (-1.356 + i 4.482) fm, 

(n) 



n = (0.1679 - i 0.6937) fm, and T^i = (0.003497 - 
i 0.0003535) fm. We conclude by noting that the d{ 



(») 



in both n- 3 He and p- 3 R arc not purely isoscalar, but and d 2 combinations give the leading contributions to 
rather have significant admixtures of isospin components J n an d that, in the case of pion exchange, d[ 1] is in 
T > 0. fact dominant. This fits in well with the expectation 
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1 oo.ii 


lp!il,L rp'21,1 

•'01,10 01,11 


n 


Re 


Im 


Re 


Im 


Re 


Im 


1 


-0.104E+01 


-0.302E+01 


-0.133E+00 


-0.134E-02 


-0.316E-01 


-0.168E-01 


2 


+0.219E+00 


-0.996E-01 


-0.830E-02 


+0.143E-03 


+0.210E-01 


+0.123E-02 


3 


-0.289E+00 


+0.108E+00 


+0.129E-01 


-0.269E-03 


+0.136E-02 


-0.293E-03 


4 


-0.767E-01 


-0.971E-01 


-0.442E-02 


+0.335E-05 


+0.401E-02 


+0.330E-03 


5 


-0.771E-01 


-0.111E+00 


-0.573E-02 


-0.682E-05 


+0.463E-02 


+0.385E-03 


6 


-0.226E-02 


-0.285E-02 


-0.625E-03 


+0.453E-04 


-0.104E-01 


-0.224E-03 


7 


-0.110E-02 


-0.210E-02 


+0.150E-03 


+0.194E-04 


-0.126E-01 


-0.277E-03 


8 


-0.393E+00 


+0.159E+00 


+0.124E-01 


-0.253E-03 


+0.189E-02 


-0.239E-03 


9 


-0.161E+00 


+0.559E-01 


+0.637E-02 


-0.144E-03 


+0.110E-01 


+0.364E-03 


10 


-0.843E-01 


-0.107E+00 


-0.495E-02 


+0.375E-05 


+0.445E-02 


+0.367E-03 


11 


-0.887E-01 


-0.126E+00 


-0.667E-02 


-0.757E-05 


+0.533E-02 


+0.443E-03 


12 


-0.265E-01 


-0.295E-01 


-0.640E-02 


-0.285E-04 


-0.340E-03 


-0.245E-03 



TABLE XIV: The parity-violating T-matrix elements corresponding to the DDH potential components O^ in combination 
with the AV18/UIX strong interaction potentials at vanishing n- 3 He energy. The statistical Monte Carlo errors are not shown, 
but are typically less than 10%. The T-matrix elements are adimensional, see text for explanation. 



n 


Redi" J 


Re4 nJ 


R»4 BJ 


Red[ n) 


Re4 n; 


rDDH 


1 


+0.617E+00 


+0.349E-01 


+0.107E-01 


-0.414E-02 


+0.616E-04 


-0.189E+00 


2 


+0.384E-01 


+0.436E-01 


+0.333E-03 


+0.345E-02 


-0.597E-05 


-0.577E-02 


3 


-0.598E- 01 


-0.559E-01 


-0.356E-03 


+0.249E-03 


+0.923E-06 


+0.864E-02 


4 


+0.205E-01 


-0.614E-02 


+0.347E-03 


+0.651E-03 


-0.147E-05 


-0.794E-02 


5 


+0.266E-01 


-0.527E-02 


+0.395E-03 


+0.751E-03 


-0.171E-05 


-0.976E-02 


6 


+0.287E-02 


-0.183E-03 


+0.102E-04 


-0.173E-02 


+0.159E-05 


-0.170E-03 


7 


-0.709E-03 


-0.378E-04 


+0.748E-05 


-0.210E-02 


+0.195E-05 


+0.115E-02 


8 


-0.573E-01 


-0.769E-01 


-0.528E-03 


+0.333E-03 


+0.694E-06 


+0.404E-02 


9 


-0.294E-01 


-0.309E-01 


-0.184E-03 


+0.181E-02 


-0.213E-05 


+0.289E-02 


10 


+0.230E-01 


-0.672E-02 


+0.383E-03 


+0.722E-03 


-0.163E-05 


-0.887E-02 


11 


+0.310E-01 


-0.612E-02 


+0.451E-03 


+0.864E-03 


-0.197E-05 


-0.113E-01 


12 


+0.297E-01 


-0.241E-02 


+0.106E-03 


-0.401E-04 


+0.888E-06 


-0.979E-02 



TABLE XV: The real parts of the coefficients di (i — 1, . . . , 5), and the coefficients /J? DH , corresponding to the DDH potential 
components O™ in combination with the AV18/UIX strong interaction potentials. The statistical Monte Carlo errors are not 
shown, but are typically less than 10%. The di are adimensional, while 7„ DH is in units of fm _1 . 



that the ^o — > 3 Po transition entering Tq ' 1:l in d\ is 
predominantly isoscalar, while d\ involves the transi 



tion 3 Si 



*1 m 1 01,10' 



which presumably has both 



isoscalar and isovector character. Indeed, the contribu- 
tions of isoscalar p- and w-exchange interactions are com- 



,(3) 



parable vn d\ , d\ and d\ , d\ , respectively. 
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